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Abstract. We present the broadband phase averaged spectrum of one of the brightest intermediate polars V1223 
Sgr, obtained with INTEGRAL and RXTE observatories (3-100 keV). Good statistical quality of the spectrum 
in a hard X-ray energy band (INTEGRAL/IBIS and RXTE/HEXTE) allowed us to disentangle contributions of 
a direct optically thin plasma emission and a reflected component to the spectrum of V1223 Sgr. The obtained 
measurement of the post-shock temperature of the accreting matter give us the information about the mass of 
the white dwarf and the inclination of the system. 
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1. Introduction 

V1223 Sgr is a very well known binary star that belongs 
to the class of intermediate polars (IPs) - binary systems 
with a white dwarf that accrets matter from a Roche lobe 
filling secondary star (see e.g. comprehensive model and 
discussion in lBeuermann et al. 2004)) . In the IP systems 
magnetic field of a white dwarf is strong enough to con- 
trol the accreting matter at relatively large distances from 
it and to form a magnetosphere and magnetic funnel (see 
e.g. lPatterson 1994l for a review). The matter falls down to 
the surface of the white dwarf with almost free fall veloc- 
ity and is heated in a strong shock to temperatures of the 
order of ~ 10 keV. Falling matter preferentially comes to 
the magnetic poles and emits X-rays via an optically thin 
thermal plasma emission (see e.g. lLamb fc Masters 19 79). 
In spite of this relatively simple scheme the observed emis- 
sion of IP is much more complicated. First of all its spec- 
trum can contain a set of components with different tem- 
peratures (mostly visible at low energies), it often demon- 
strates an intrinsic photoabsorption, partial covering, dif- 
ferent ionization states of an emitting plasma and a re- 
flection from an optically thick relatively medium (see 
e.g. I Warner 1 995 for a review). Often the spectrum de- 
pends on orbital phase and phase of rotation of the white 
dwarf (e.s. INorton fc Watson 1 989). Different objects of 
this class were extensively studied by several X-ray obser- 
vatories over last decades. However, until now there was 
no study of hard X-ray emission of IP (E > 20 keV) except 
for the case of V709 Cas (de Martino et al. 2001). 
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Energy band > 20 keV is important because it allows 
to confidently measure the maximal temperature of the 
shock in the accretion column of the IP. It is especially 
important if the X-ray spectrum of IP contains strong re- 
flected component, that distorts the observed continuum 
at energies ~10-20 keV. In the most cases the temper- 
ature of the post-shock region of the accretion column 
of IP is of the order of ~ 10 — 40 keV that makes it 
quite difficult to measure with telescopes working at 0.5- 
10 keV energy band, like EINSTEIN, ASCA, CHANDRA, 
XMM. Observations of IP performed with GINGA/LAC 
(e.g. Ishida 1991, Ishida et al. 1994) practically do not 
provide statistically significant measurements at E ~20 
keV, that makes it hard to confidently measure the tem- 
peratures higher than ~ 10 — 20 keV. As the energetics of 
the X-ray spectrum of IP is dominated by bremsstrahlung 
emission with maximal temperatures the accuracy of its 
measurement is of high importance (see. e.g. discussion in 
IBeuermann et al. 2004(1 . 

In this work we present results of measurements of the 
broadband spectrum of the intermediate polar V1223 Sgr 
with the RXTE and INTEGRAL observatories with an 
accent on the averaged hard X-ray part of the spectrum. 



2. Data analysis 

The international gamma-ray observatory INTEGRAL 
was launched by the Russian launcher PROTON from 
the Baikonur cosmodrome in the high-apogee orbit on 
October 17, 2002 (Eismont et al. 2003). The payload in- 
cludes four principal instruments which allow one to carry 
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Fig. 1. X-ray map of the region around VI 223 Sgr ob- 
tained with INTEGRAL/IBIS /ISGRI in the energy band 
18-60 keV 



out simultaneous observations of sources in the X-ray, 
gamma-ray and optical energy range (Winkler et al. 2003). 

In the August-September of 2003 ultra deep (~2 mil- 
lion sec) observations of the Galactic center region were 
performed. The survey of this region with catalog of de- 
tected sources was published by Revnivtsev et al. (2004). 
The intermediate polar VI 223 Sgr was detected during 
these observations with a high statistical significance in 
the 18-60 keV energy band. 

Due to relatively large distance from the Galactic 
Center (~ 16°), V1223 Sgr was outside the field of view 
of the X-ray monitor JEM-X of the INTEGRAL obser- 
vatory. On the other hand the source is weak and rela- 
tively soft, therefore the sensitivity of the spectrometer 
SPI l |Roques et al. 20031 ) is n °t enough to measure the 
spectrum of the source. Therefore in the paper we will 
concentrate on the data obtained by IBIS telescope, its up- 
per detector ISGRI IjLebrun et al. 2003|) . The source was 
quite far from the center of the field of view of the IBIS 
telescope, resulting in a strong (factor of ~9.5) reduction 
of its effective area useful for the measurement of the spec- 
trum of V1223 Sgr. 

Data of the IBIS telescope were processed using meth- 
ods described in Revnivtsev et al. (2004). The image of 
the region around V1223 Sgr is presented in Fig.l. For the 
spectral analysis we used a simple ratio of the flux mea- 
sured from the source to the fluxes of the Crab nebula in 
the same energy bands and assumed that the Crab nebula 
spectrum has the form dN(E) = \QE~ %l dE, where N(E) 
is the number of photons at energy E. As the spectrum 
of V1223 Sgr does not demonstrate very sharp structures 
such method of the spectral reconstruction is acceptable. 
Tests of this method on observations of the Crab nebula 
showed that the resulted spectra can still possess approxi- 



mately 5% systematic uncertainties in every energy chan- 
nel. 

In order to compare the obtained results with 
measurements of other instruments and to accomplish 
the obtained hard X-ray spectrum of the source with 
standard X-ray band measurements (~l-20 keV) we 
used data of the RXTE observatory (jBradt et al. 1 993 
IRothschild et al. 1998j> . The PCA and HEXTE spectrom- 
eters of this observatory cover the energy range 3-250 
keV. Effective areas of these spectrometers are quite large 
~ 6400 cm 2 for the PCA at energies -6-7 keV and ~650 
cm 2 for each (of two) clusters of HEXTE detectors at en- 
ergies ~50 keV. We quote here effective area of only one 
HEXTE cluster of detectors because only one HEXTE 
cluster can observe the source at the same time, the other 
one observes the background at this time. 

Intermediate polar V1223 Sgr was observed by the 
RXTE observatory 14 times during period 1996-2000. 
Total usable exposure time of these observation is approx- 
imately 114 ksec. 

The data of the RXTE observatory was reduced with 
the help of standard LHEASOFT/FTOOLS 5.3 package. 
As the spectrum of IPs does not vary much at ener- 
gies higher than >3 keV (see e.g. iBcardmo re et al. 2 000 
for the case of V1223 Sgr) for the following analysis 
we averaged all available data and analyzed only aver- 
aged spectrum of V1223 Sgr. For the spectral approxi- 
mation we used XSPEC package. Normalizations of spec- 
tra obtained by different instruments were set free to 
vary during the fitting procedure. Systematic errors 1% 
and 5% were added to every channel of the RXTE/PCA 
and INTEGRAL/IBIS spectra, respectively. The resulted 
spectra of all instruments were rescaled to that of the 
RXTE/PCA. 

3. Spectral analysis 

First of all we approximated the hard X-ray part of the 
spectrum of V1223 Sgr obtained with INTEGRAL/IBIS 
and RXTE /HEXTE (energy band 20-100 keV) with 
a simplest thermal model - an optically thin thermal 
bremsstrahlung. The quality of the fit is very good (re- 
duced x 2 ~ 0-9) and the parameter of the temperature 
is equal to kT = 29 ± 2 keV. Best fit parameters for 
separate data of RXTE/HEXTE and INTEGRAL/IBIS 
are 25.6 ± 2.5 keV and 29 ± 5 keV, respectively. As the 
obtained results agree well with each other in our subse- 
quent analysis we will combine data of these two instru- 
ments. Obtained value of the temperature parameter is 
somewhat smaller than that obtained by Beardmore et al. 
(2000) from GINGA data (kT = 43 ± 13 keV). However 
if we apply the same model (with the neutral absorption 
and broad gaussian line in the energy range ~ 6 — 7 keV) 
to the data of the RXTE/PCA detector (energy band 3- 
20 keV) not taking into account harda X-ray data the 
resulted temperature would be in a very good agreement 
with that of Beardmore et al. (2000) - kT = 44.5 ± 1.5 
keV. Moreover, if we will try to approximate only spectral 
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Table 1. Best fit parameters for averaged broadband (3- 
100 keV) spectrum of V1223 Sgr 



V1223 Sgr 



Parameter 




N H , cm" 2 




3.3 ±0.2 


kT, keV 




29 ±2 


R = Q/2n 




0.35 ±0.15 


-E^gauss 




6.63 ±0.03 


& gauss 




0.37 ±0.02 


EW, eV 




387 ± 18 


FluX obs , 10" 10 


erg/s/cm 2 


3.76 ± 0.04 a 


Flux corr , 1(T 10 


erg/s/cm 2 


4.57±0.05\ 


X 2 /d.o.f. 




76.5/89 



- R — Q,/2tx - is the parameter of the reflection, solid angle 
subtended by a reflector 

a - observed flux in the energy band 3-100 keV, the flux was 
calculated assuming that the spectrum of the Crab nebula 
has the following normalization dN(E) = 10.0£~ 21 <±E. If 
the assumption of the Crab nebula spectrum is different 
the flux value should be rescaled 

b - absorption corrected flux in the energy band 0.5-100 
keV 

points in the energy band 10-20 keV we will obtain even 
higher temperature - kT = 65 ± 5. This discrepancy can 
be understood if we will take into account that the spec- 
trum of V1223 Sgr contains relatively strong reflected (e.g. 
Bas ko, Sunyaev fc Titarchuk 1974| ) component and there- 
fore the model for the spectral approximation should be 
more complex. 

In the work of Beardmore et al. (2000) only 
GING A/LAC (energy range 2-37 keV) data were used 
to describe the source spectrum at energies >10 keV. 
However, the spectrum of V1223 Sgr obtained by 
GING A/LAC practically does not contain statistically sig- 
nificant points at energies higher than ~ 20 keV . In 
spite of the fact that Beardmore et al. (2000) included 
reflected component into their spectral approximation it 
is very hard for their dataset to determine at the same 
time the temperature of the thermal plasma and the am- 
plitude of the reflection. In our case we have data of 
INTEGRAL/IBIS and the RXTE/HEXTE (20-100 keV) 
and can determine this parameters more confidently while 
fitting the broadband source spectrum. 

Approximation of the broadband spectrum 
(BeppoSAX data) of the another intermediate polar 
V709 Cas made by de Martino et al. (2001) confirms 
these arguments. Authors showed that the broadband 
spectrum of IP should be described by a thermal model 
with the reflection. Neglecting the X-ray reflection in 
a spectral modeling result in the overestimation of the 
temperature of a post-shock plasma in the accretion 
column of IP. 

For the broadband spectral approximation we used a 
simple analytic model consisting of the bremsstrahlung 
emission with its reflection from an optically thick cold 
medium plus the neutral absorption and broad gaussian 
line at ~6.5 keV. The latter component mimics a set of 
lines of optically thin plasma and fluorescent Fe line, orig- 
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Fig. 2. Averaged spectrum of V1223 Sgr obtained with 
INTEGRAL and RXTE observatories. Solid curve repre- 
sents the best fit model. Dashed line shows the emission 
of the thermal plasma, dotted line - the contribution of 
the reflected component. On the lower panel the differ- 
ence between data and the model in the units of standard 
deviations (a) is shown. 
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Fig. 3. Confidence contours (1,2 and 3 a) of the plasma 
temperature (kT) and relative reflection amplitude (R = 
Q/2tt) 



mating as a result of the reflection. Such spectral model 
should approximately represent emission of the hot plasma 
in the post-shock region of V1223 Sgr. The best fit tem- 
perature of the bremsstrahlung emission in this case will 
denote the maximal temperature in the post-shock region 
of the accretion column of IP. Obtained best fit param- 
eters of the approximation of the combined INTEGRAL 
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and RXTE spectrum of VI 223 Sgr are presented in Table 

1. The averaged spectrum of V1223 Sgr is shown in Fig. 

2. As the value of the plasma temperature strongly corre- 
lates with the amplitude of the reflection we also present 
confidence contours on these values in Fig. 3. 

4. Discussion 

In the previous section we analyzed the broadband phase 
averaged spectrum of one of the brightest intermediate 
polar V1223 Sgr. We have shown that even apart from 
detection of fluorescent Fe line at energy ~ 6.4 keV (see 
IBeardmore et al. 2000|l there are strong indications on the 
presence of the reflection component in the spectrum of 
V1223 Sgr. Because of this component any measurements 
of the post-shock temperature of the accretion column 
near the surface of IP that are based on instruments with 
an effective energy band lower than ~20-30 keV are bi- 
ased towards obtaining higher temperatures. Our spectral 
modeling for the first time includes hard X-ray energy 
range (20-100 keV) and provides more confident estima- 
tion of the post-shock temperature. Using this value we 
derive some important parameters of V1223 Sgr. 

Current models of a shock in the accretion column of 
IP assume that the accreted matter is heated up to its 
virial temperature (see e.g. lAizu 19731 IFabian et al. 19761 
ILamb & Masters 1979|l . The heated plasma produces 
mainly bremsstrahlung with possible contribution from 
a cyclotron emission at low accretion rates (see e.g. 
IFabian et al. 19761 IWarner 1995jl . 

Our new broadband measurements of the X-ray spec- 
trum of VI 223 Sgr provide us a set of observables with 
higher confidence than that of previous. They are: 1) the 
post-shock temperature T s = 29 ± 2 keV, 2) the broad- 
band X-ray flux (absorption corrected) Fo.S-iookeV = 
4.6x 10 -10 erg s _1 cm -2 , that in turn gives us the emission 
measure of one emitting accretion column i? m .i = J N^dV 
= 1.07 x 10 13 47T.D 2 cm -3 , where N e - is the number den- 
sity of electrons in a hot emitting plasma. 

The distance to V1223 Sgr is now known from HST 
measurement of the parallax of the system, D = 527^|g 
pc (Beuermann et al. 2004). Therefore we can evaluate 
the isotropic luminosity of one accretion column L x .i = 
47rL> 2 J F u . 5 _i 00kcV = (1.5 ± 0.3) x 10 34 erg s"\ and the 
volume emission measure E mi x ~ 3.4 x 10 56 cm -3 . 

The shock temperature can give the direct evalua- 
tion of a white dwarf mass (see e.g. IFabian et al. 19761 
IWarner 199"5|> 

r , = 3 l.B^bV. 

Mq i? W D 

Using an approximation of the white dwarf mass- 
radius relation from Nauenberg (1972) 



Ro = 0.78 



1.44M 



M 



2/3 



M 



1.44M, 



2/3' 



1/2 



where Rg - is the radius of the white dwarf in the units 
of 10 9 cm, we obtain for the mass and radius of the white 
dwarf in the V1223 Sgr: 

M WD = 0.71 ± O.O3M 

R WD = (7.8 ± 0.2) x 10 8 cm 

The obtained estimation of the white dwarf mass is 
slightly smaller than that obtained by Beardmore et al. 
(2000) because of lower post-shock temperature measured 
by us and is better consistent with the evaluation of 
white dwarf mass from optical observations: A/wd = 
0.4 - 0.6M© dPenning 19851 

The mass function of the secondary star in the system 
V1223 Sgr is f(M 2 ) = 0.0026 ± 0.0015 M Q (Watts et 
al. 1985), the mass of the secondary star is M 2 = 0.4 
Mq (Penning 1985). Using these values we can obtain the 
inclination angle of the system i = 22° ±5°. If we will 
assume that the mass of the secondary star is lower, e.g. 
Mi = 0.25M© (|Beuermann et al. 2004j> the estimation of 
the inclination angle will increase to i = 32° ± 7° . 

As the X-ray emission is the major channel of energy 
losses in the post-shock region we can derive the mass 
accretion rate M 



M 



3 /j,mn2L x , 
8 kfi 



where /i=0.62 is the mean molecular weight of the accret- 
ing matter. Here factor 2 stands in order to account for 
two emitting magnetic poles of IP. Substituting measured 
values we obtain M ~ 2.54 x 10 17 g s -1 . 

Using M WD = 0.71M Q , R wd = 7.8 x 10 s cm and 



E mA = 4tt.R wd //2 h s Nj ~ 3.4 x 10 56 cm" 3 , where / - is 
the fraction of the white dwarf surface, occupied by two 
accretion columns, we can evaluate the height of standing 
shock h s and the number density of electrons in the post- 
shock region N e following Fabian et al. (1976) and Warner 
(1995). 



N e = 3.1 x 10 13 M 17 



3.3 x 10 4 / Myj 



WD 



Mr. 



1/2 



R 



-3/2 J.-1 



( M- 



WD 



3/2 



R^ 2 cm 



where Myj - is the mass accretion rate in units of 10 17 
g s -1 . We can derive N e — 9.6 x lO 16 /^ cm~ 3 and 
h s = 1.1 x 10 7 /_3 cm, where /_3 = 10~ 3 / - is the fraction 
of the surface of the white dwarf, occupied by two accre- 
tion columns in the units of 10 -3 . Note that the derived 
parameters are similar to those obtained by de Martino et 
al. (2001) for the intermediate polar V709 Cas. 
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